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1 Introduction 

Motivation. Probabilists share with physicists the tradition of paying special 
attention to analytically tractable models, like for example solvable Markovian 
semigroups/evolution equations, in their respective parallel terminologies. Such 
models have lead to significant progress over the centuries, as witnessed for exam- 
ple by the recent cases of the Black-Scholes and affine models, in mathematical 
finance. 

With a few exceptions, the main approach used nowadays for dealing with 
Kolmogorov evolution equations (to adopt a 'unified' name) in applied proba- 
bility is the transform method initiated by Laplace, Poincare and Weierstrass. 
It is intriguing to ask whether the alternative algebraic-geometric approach pi- 
oneered by Lie [15], which is going these days through a revival period (see, for 
instance, |8I10I12I18I25| ) stimulated by today's symbolic computing, could turn 
out to be as useful in probability. 

Note that Lie's approach would seem to be taylor made for jump-diffusion 
processes, which are generated by combining noncommuting operators. 
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Jump-diffusions. A recurring theme in applied probability is distinguishing 
between two possible sources of uncertainty: small continuous changes modeled 
by diffusion, and 'catastrophic' changes modeled by a jump process. To resolve 
this issue, one uses jump-diffusions models, i.e. solutions of a SDE (stochastic 
differential equation) 

dX t = <p{X t )dt + o{X t )dB t - dS t , (I) 

where 

— B t is standard Brownian motion, 

— St is a pure jump process, with a Levy density v(x, z) := X(x) b{z), which may 
arise, for example, from i.i.d. jumps Cj (sometimes one sided, for example 
negative) , whose density, distribution, complementary distribution, and first 
moment are denoted respectively by b(x), B(x), B{x), b\. 

The first two terms (the drift ip and the variance er) of the Levy-Khinchine 
triple ip, er, v define a continuous diffusion process, and the last term defines a 
pure jump/convolution process. 

Jump-diffusions (JTJ) are Mar kovian processes with associated evolution/ back- 
ward Kolmogorov equation 

WM = g xfM , f(x,0) = f Q (x), (2) 
which describes expectations evolving in time 

f(t,x)=E Xo=x f Q (X t ). 
The infinitesimal generator operator is given by 

Qf(x) = G«f{x) = <p(x)f'(x) + ^P-f"(x)+ 

{f{x -z)- f{x))v{x, z)dzQ^f(x) + g®f(x), 



for any twice continuously diffcrentiable and bounded function f(x), where the 
second part Qx^ is associated to the pure jump convolution part. 

One important example, already studied in Kolmogorov's founding paper 
[IS] , is that of hypergeometric diffusions with quadratic variance and linear drift: 

d 2 f df 
G {d) f(x) := (a 2 x 2 + a\x + a ) ^-j + ((fix + ip ) q^- 

The Levy model is obtained when the variance and drift rates as well as the 
Levy intensity v (x, z) are independent of x. 
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Example 1. The Cramer Lundberg risk model (1903) is one of the most studied 
models in applied probability |20j , describes the surplus of an insurance company: 



with initial capita u, linear premium rate/drift ct and claims Ck, modeled by a 
sequence of i.i.d. positive random variables with a common density f(x) = fc{x), 
and which arrive at the increase points N — {N t ,t > 0} of an independent count- 
ing process with EJVj = Xt. If moreover Nt is a Poisson process with exponential 
interarrival times, then S(t) is a compound Poisson process with positive sum- 
mands, and U(t) 1 with t € R+, is Markovian. 

If, more generally, S(t) is a subordinator (nondecreasing Levy process), then 
U(t) is a spectrally negative Levy process (with jump measure supported on the 
negative half- line). 

Phase-type distributions. For a general line of attack on jump diffusions, 
it seems natural to start by adding an independent compound Poisson process 
with phase- type jumps C, (the convolution term in the evolution equation), with 
distribution 



where B,/3, 1 are respectively a subgenerating matrix (nonnegative off-diagonal 
elements and nonpositive row sums), a probability row vector and a column 
vector of ones. Note that the density (the negative of the derivative of B(x)) is 



where b = — B 1 is a column vector with n components. 

Phase- type distributions have the advantage that the intcgro-diffcrcntial equa- 
tions satisfied for example by the ruin probabilities ^(t, x) may be transformed 
into ordinary differential equations -see for example Paulsen [22] , which resolves 
by this approach four notoriously difficult particular cases, in terms of special 
hypergeometric functions. Paulsen's paper arose the natural question of whether 
under more general models of hypergeometric diffusions with phase-type jumps, 
ruin probabilities could be determined analytically as well. 

Some further particular examples were studied for example in the papers 
[4 13]. We are continuing here this line of research, providing a new solvable 
model. 

Note that it is also possible to consider the more general matrix-exponential 
distributions, where one assumes only the nonnegativity and integrability (to 1) 
of the density b(x), but the advantages are not clear, since phase-type distribu- 
tions are already dense [3]. 

First passage problems. Denote by 



N(t) 




B(x) := P[C >x] = /3e Bx l, 



b(x) = /3e Bx b, 



.+ 

L 



mf{t > 0; X t > L} 



t = n= inf{f > 0; X t < 1} 
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the first passage times of a stochastic process above/below given levels L, I. 
The latter, also called ruin time in the insurance literature, is one of the old- 
est applications of probability and ODE's, introduced by Thiele, the founder 
of the Danish insurance company Hafnia (1872) -see www.stats.ox.ac.uk/ stef- 
fen / seminars / centertalk.pdf. 

Ruin probabilities. The first objects of interest in first passage theory are 
the finite-time and ultimate/infinite horizon ruin probabilities W(t,x) and the 
related survival probabilities \P(t, x) 

&{t,x)=P x [T<t], &{x) = P x [t < oo], 

W(t,x) =P x [t >t] = l-W(t,x), W{x) =P x [t = oo]. 

For the Markovian first step/infinitesimal analysis shows that the 

ultimate ruin probabilities are harmonic functions, satisfying 

gW{u) := — W"(u) + c&'(u) - M> (u) + A / &{u- z)f(z)dz + XF(u) = 0, (3) 
2 Jo 

with fy(u) — 1 and u < 0. 

Note 1. Note that the trivial solution of equation ([3]) with fy(u) = 1 (true when 
lim^oo X{t) = — oo) may be discarded by adding the restriction 

lim X(t) =ooo lim &(u) = 0. 

t— >oo u— >oo 

The condition for that to hold is the same as for the process reflected at 
to be nonergodic; in the Levy case, this is quite simple: 

E ^)-x(o )=c _ AECi>0) 

(by the classic law of large numbers), but things get more complicated for the 
generalizations considered here, with the nonergodicity condition involving the 
invariant measure P(x) of the process, i.e. the nonnegative solution of the adjoint 
equation Q* p(x) = 0. 

Note 2. In ([3]) and in any problem involving the process under consideration 
must appear the same operator Q, which is also the generator of the associated 
semigroup of transition operators. Put informally, Q is the 'analysis dictionary 
entry' associated to a given Markovian process. For example, the finite-time 
ruin probabilities must satisfy the nonautonomous backward Kolmogorov partial 
intcgro-diffcrential equation: 

d 

- —&(u,t) + g&(u,t) = 0, &(u, t) = 1, u < 0, Vt. (4) 
at 

Laplace transform in time. The finite-time ruin probabilities are very 
seldom available analytically; one natural way to approximate them is by com- 
puting and inverting their Laplace-Carson transform 

I* oo 

V q {x) = / qe- qt V{t,x)dt 7 (5) 
Jo 
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which may be also viewed as the probability of passage of the process 'killed' 
after an independent exponential random variable e q of rate q > 0: 

¥ q (x) = Pr(r < e q ) = E x { e -" T I {t<oo} ) = 

/ e _ «* Pr(r £ dt) = / qe~ qt 'P(t,x)dt, 
Jo x Jo 

where the last equality follows by integration by parts. 

More generally, we want to calculate the killed ruin probability 

$ q (x)=* l q (x)=E x e-fo*W*))<i°, ( 6 ) 
for general discount functions q(X(s)), and the killed ruin with overshoot penalty 

= E a (e-/o T <KX«)^(* T -<) /{T<oo}) . 

By Dynkin's formula, this function is a solution of the Sturm- Liouville type 
equation 

Gf(x)-q(x)f(x)=0, (7) 

f[x) =e^- l \x<l, f(l) = l, ifa^O, (8) 

see for example [6j pag. 824]. Note that the number of boundary conditions 
depends on the presence of Brownian motion (probabilistically, this is required 
by the presence of a new unknown: the probability of crossing continuously). 

Note 3. The arguments I, £ will be often suppressed (in particular when they 
equal 0). 

The evolution equation ([2]) and the corresponding time-independent coun- 
terparts, the invariant measure and the harmonic functions of interest in first 
passage theory, have been intensively studied for diffusions and for Levy pro- 
cesses. 

In both the diffusion and the Levy case there exists a beautiful first-passage 
theory. In the latter case, this reduces the calculation of many of the functions 
(0 to the Wiener-Hopf factorization (of the generator, or rather of its symbol) 
-see for example |5I16| . 

A question which begs of itself is combining these two well understood gener- 
ators, i.e. investigating the existence of extensions under some common umbrella 
case like processes with rational generators. Note that the affine case has already 
been thoroughly investigated, for example in mathematical finance, for modeling 
interest rates. 

Below, we focus on the case of piecewise deterministic processes (no diffu- 
sion) with phase-type downward jumps. In this case, the Feynman-Kac integro- 
differential equation |14] for killed ruin probabilities may be brought to the form 
of a ODE linear system [3] of the form 

/ &'(x)\ /a+£ z$0\ , , 
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where B is a n x n stochastic generating matrix, where M is a column vec- 
tors with n components, and (3 = (J3i, . . . ,/3„) is a row probability vector, i.e. 
S"=i Pi = 1- The variable W is the killed ruin probability, the function (p(x) is 
the drift, the constant q is the killing rate/Laplace transform argument, and the 
components M\, . . . , M n of the vector M are killed ruin probabilities, obtained 
by changing the jumps to segments of slope ±1 for upwards/downwards jumps, 
and by associating 'auxiliary stages of artificial time' to the phases of the jumps 
(Mi(x) is thus the killed ruin probability when starting at x in phase i). 

Note that with upward jumps, the analog equation will have the last rows 
corresponding to the M variables multiplied by —1. 

Note 4- We will restrict to first-passage problems in domains where the drift 
f(x) does not change sign, which determines then corresponding boundary con- 
ditions. 

The theory of Lie systems |8I10|19I25] , of which nonautonomous systems of 
first-order linear homogeneous differential equations are a particular case, states 
that linear systems like © are integrable by quadratures, i.e. their solutions 
can be explicitly obtained in terms of 'algebraic operations' and integrations of 
certain given functions depending on one variable |17j . if there exists a finitely 
generated solvable matrix Lie algebra g such that 

1 = f ^ I =-rrf + ^) efl, V^i. (10) 




<p(x) V ) \h B 

For example, in the easy case of a constant drift <p(x) = c, the integrability 
condition (flOl) is trivially satisfied as A x 6 g = (Aq), for all x £ K, and the 
solutions are phase-type functions. 

In the case of a non-constant drift <p(x), verifying condition (fTU)) reduces to 
proving that the sequence of commutators 



[Tx,T 2 ], [Tx, [Tx,T 2 ], [T 2 , [Ti, T 2 ]], [T x , [. . . [T^] . . .}],. . . (11) 
generated by the matrices 




Tx = , f 2 = 




b B 

must contain a finite subset (generators) spanning a finite-dimensional Lie alge- 
bra with respect to the matrix commutator. 

The example we focus on below is that of downward exponential jumps of 
rate \x, over an exponential horizon e q , when the linear system ([9]) becomes: 

M'(x))-{ y -y)\M(x)) (12) 

with /i>0, A>0, q > (here (3 = 1, the unique 'probability vector' of 
dimension 1). 
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Example 2. When q — 0, the solution to system (fT2|) satisfying the conditions 
<^(oo) = Af(oo) = and M(0) = 1 can be derived by subtracting its equations, 
yielding: 

&(x) - M(x) = (if (0) - l)e z( - x \ Z(x) = -fix + [ -^—dv, 



o <f(v) 

and, provided Z{od) = — oo, then 

M(x) =n(l-V(0))J e z{v) dv^^{x) = (1-^(0)) (u J e z ^dv - e z{x) ^j . 

With <p(x) = c constant, we have M{x) = e~( 1 ~ v ') ,ix , where r\ is the smallest 
positive root of cr)(—fi + fit]) — ((A + q)rj — A) = 0, and when q — 0, we have 
rji = ^j, i]2 = 1) from where we recover the well known &(x) = *~ • 

We will see below that when q = 0, the theory of Lie system shows that the 
above system satisfies the integrability condition (fTOj) . and it can be therefore 
integrated by quadratures. Otherwise, Theorem [5] ensures that this is no longer 
the case. 

It is natural to ask whether other families of processes with solvable first 
passage probabilities exist. By the alternative Riccati approach carried out in 
Section |H we take advantage of the fact that several families of Riccati equations 
integrable by quadratures have been accumulated in the literature, like the gen- 
eralised Allen-Stein family described in Theorem [TJ Hovy's equation obtained 
with ip{x) — x,q — (in fact, our equation is a 'nonlinear Hovy's equation'), or 
Robin's family [33]. Each of these cases will lead us immediately to a class of 
piecewise deterministic processes with exponential jumps which can be solved 
analytically. We exemplify this for the generalised Allen-Stein family, formed by 
those particular cases of the system (fL2|) satisfying the integrability condition 



(p'/2 + (X + q)- fitp = kci v /-/iAc c 2 </J, 

for certain reals constants Co, c\, c%, and with k = sg(A/i/ (ip(x)coC2)) ■ 

Note that the unifying general approach based on the theory of Lie systems 
recently provided in [10] should allow recognizing and solving higher dimensional 
first passage problems. 

Contents: Our paper illustrates how the theory of Lie systems allows recov- 
ering known results and provide new ones for first passage problems of jump- 
diffusions with phase-type jumps. 

We consider various particular instances of piecewise deterministic processes 
with exponential jumps, both associated to solvable and, more interesting, to 
non-solvable Vessiot-Guldberg Lie algebras in Sections [3J and [5] 



2 Fundamentals of Lie systems 

Let us briefly recall the most fundamental notions of the theory of Lie systems 
and the geometrical treatment of differential equations to be used throughout 
our work. 
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A fundamental concept in the geometric study of differential equations and 
in what follows is the hereby called x-dependent vector field. This concept refers 
to a map 

X : R x N — > TN 

(x,y) ^ X(x,y)eT x N, 

such that X x : y E N H> X x (y) = X(x,y) G TN is a standard vector field on N, 
for every Note that, in consequence, giving a x-dependent vector field is 

equivalent to providing a family of vector fields {X x } xe s. on N parametrized by 
what explains its name. 
In similarity to standard vector fields, each x-dependent vector field, e.g. 
X(x,y) = Yl7=i ^ l ( x > y)d/dx z , also admits integral curves (see [7]) determined 
by the solutions of the system 

^£=X\x,y), i=l,...,n. (13) 

Indeed, the relevance of studying the properties of x-dependent vector fields 
mostly relies on investigating systems of the above form by means of the analysis 
of the properties of its corresponding X(x,y), its denominated associated vector 
field. 

The theory of Lie systems mainly deals with a class of systems of first-order 
differential equations, the so-called Lie systems, associated with x-dependent 
vector fields of the form 

r 

X(x,y) = Y / b a (x)X a (y), (14) 

a=l 

where X\ , . . . , X r are a set of vector fields on N spanning a finite-dimensional 
Lie algebra of vector fields V, the associated Vessiot-Guldberg Lie algebra. These 
vector fields can be, at least locally, considered as the fundamental vector fields 
of a certain action <L> : G x N — > N, where G is a Lie group with Lie algebra 
isomorphic to V. In those cases where the form of ^ can be explicitly determined, 
the general solution, y(x), of system (fl~3f can be cast into the form y(x) = 
4>(g(x),yo), with j/o being any point of N and g{x) being the solution of the 
equation on G given by 

^- = -J2b a (x)X*(g), g(0) = e, (15) 

a—l 

where the X^ are certain right-invariant vector fields on G, for details see [8118] . 

A particular instance of Lie system having a special relevance to our work is 
given by the following system of homogeneous linear differential equations 
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associated with the X- dependent vector field X(x,y) = Yll j=i A % j{x")y^ d J dy 1 . 
In fact, consider the family of vector fields of the form {Vy — y^d/dy 1 \ i,j = 
closing on the commutation relations [Xi^Xti] = S\X^ — fiiXn, for 
i,j,k,l = l...,n, see [S]. The vector fields X^ span a Lie algebra of vector fields 
V, and taking into account that X x = Y17j=i Aj(x)Xij, ^ follows that (fTo) is a 
Lie system related to the Vessiot-Guldberg Lie algebra V. 

It is important to note that some linear homogeneous systems can be con- 
sidered, additionally, as a Lie system related to another Vessiot-Guldberg Lie 
algebra Vq CC V. In order to illustrate this claim, let us consider the linear 
homogeneous system 

dy 1 Ai,^..x , a\,^.2 



dx 



^(x)y 1 + A\{x)y\ 



describing the integral curves of the x-dependent vector field 

X. = A\ { x) U° - f°) + Al(x)y>° + Afatf* 



dy 1 dy 2 J dy 1 dy 

In this case, taking as Vq the Lie algebra 

T ^ / 1 d 2 d 2 9 1 9 

v o = (y tt ~ y t~2 ' y Tn>y «~2 

\ oy L ay 1 oy L oy z 

it turns out that X x € Vq CC V for all x £ R. Consequently, system p7|) is also 
a Lie system related to the Vessiot-Guldberg Lie algebra Vq. 

There is an alternative way to determine the Vessiot-Guldberg Lie algebras 
associated with systems of the form fp~6|) . In order to do so, let us cast these 
systems into the matrix form 

A\{x) ...A\{x) 







(!) 


-I 


^ dx ' 






and consider the map p : g[(n, M) — > V satisfying that p(Mij) = — X^, with 
(Mij)f = Sftf and i,j = 1, . . . ,n. As [My, M kl ] = <S£M U - <5.jM kj , for i,j, k, I = 
1 . . . , n, it follows that the map p is a Lie algebra isomorphism. In consequence, 
if X x £ Vq for every i£8, then A(x) belongs to the matrix Lie algebra p~ 1 (Vo). 
Similarly, if A(x) € g for every x € M, then X x is contained in the Lie algebra 
p(g). Trivially, the next proposition follows. 

Proposition 1. System jib]) is related to a (solvable) Vessiot-Guldberg Lie al- 
gebra Vq if and only if there exists a (solvable) matrix Lie algebra g C g[(n,M) 
such that A(x) e g, for every iff. 
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The theory of Lie systems establishes that those Lie systems related to solv- 
able Vessiot-Guldberg Lie algebras, and whose associated action <P is expressed 
in terms of elementary functions, are integrable by quadratures. Otherwise, our 
methods, and in general all the other ones found in the literature, requires the 
Lie system to hold some kind of extra condition to be still integrable by quadra- 
tures. This fact will be used in next sections in order to integrate certain difussion 
processes (TT2|) which cannot by related to solvable Vessiot-Guldberg Lie algebras. 



3 Lie systems and integrability of piecewise deterministic 
processes 

Our goal in this section is to study the integrability properties of piecewise 
deterministic models of the form (j 12[) . In particular, our aim is to show that the 
key cases that have been exactly solved in the literature for any drift f{x) (those 
with q — studied by Segerdahl and Paulsen [22124] ). satisfy an integrability 
condition described by the theory of Lie systems ensuring that, actually, their 
solutions can be explicitly written down. Finally, our second objective is to show 
that if q ^ 0, the previous integrability condition does not hold and each case 
determined by a drift ip(x) must be analysed separately. 

The theory of Lie systems states that integrating a Lie system, e.g. one of 
the form (fl3|) . can be, at least locally, reduced to solving an equation of the form 
p3|) on a Lie group G, provided the action related to the Lie system is known. 
Moreover, it is also know [5] that if some Vessiot-Guldberg Lie algebra related 
to the Lie system is solvable, the equation (fl5|) can be exactly solved by means 
of an algorithmic method and, through the action <£, the general solution to the 
initial problem can be recovered [3]. 

Taking into account the above considerations, it is interesting to study when 
the Lie system (TT21 is associated with a solvable Vessiot-Guldberg Lie algebra 
and their solutions can be, consequently, obtained. In view of Proposition!!] this 
fact depends on the existence of a solvable matrix Lie algebra q including the 
matrices 

A ^{f '-t)=^(t"o 1 ) + "(?-i)' 

Note that when q = is assumed, the above family reads 
where 



J ' z 

The above matrices satisfy [Xi,T2] = —T\ — T 2 and span a two dimensional 
solvable Lie algebra V = (Ti,T 2 ). In consequence, system (TT2")) . with q = 0, can 
be easily solved for every drift f(x). In this way, it is not surprising that the 
solutions for these models are known since a long time ago [22|24j . 
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Proposition 2. When q =/= 0, and for a non-constant drift <f(x), the matrices 
A x span the non-solvable Lie algebra gl(2, M) of 2 x 2 real matrices. 

Proof. It is obvious that A x — X/tp(x)Ui + IJ.U2, where 

A+2 _1 ' 



In terms of the above matrices, we define 

/ 1 \ /() _> N 

C/ 3 = ([C7 1) C/ 2 ]+C/ 2 + U 1 )X/q + U 2 =\ , ?7 4 = [?7i, Z7 3 ] = 

Consequently, for every Lie algebra g such that {Ajjzgn C g, the matrices 
Ui, 1/2,1/3,1/4 must be contained in g, as they are made up from Lie brackets 
and linear combinations of elements of g. Moreover, as q ^ 0, the matrices 
Ui, 1/2,1/3 and U4 are linearly independent and they span g[(2,IR). It follows 
that g[(2,R) C g. Consequently, the Lie algebra g is not solvable. 

The consequence of the above proposition is clear: as the Lie algebra g is not 
solvable, there exists no general method to solve system (|12[) for an arbitrary 
drift when q 7^ 0. Nevertheless, despite the absence of a general method for 
solving such systems, the theory of Lie systems also provides criteria to ensure 
the integration of certain of these systems, generally satisfying some kind of extra 
condition. The applications of one of these criteria will be the main purpose of 
the following sections. 

It is important to remark that, although there are other alternative meth- 
ods to obtain the results described within this work, the theory of Lie systems 
provides a unifying approach to derive and analyse them. Moreover, we think 
that this theory gives us an appropriate approach to the analysis of systems 
(JHJ) , as it can provide generalizations of the methods here described to deal with 
higher-dimensional cases. 



4 The Riccati approach 

In order to simplify our treatment of systems (1121) , we here accomplish an alter- 
native approach to study (|T2"j) consisting on writing the system in the coordinate 
system {rj = 'P/AI, M}, bringing it to the form 

V da; 

The above non-linear system is made up from a homogeneous equation in the 
variable M and a Riccati equation in the variable 77 (with no dependence on the 
variable M), which will be called below Segerdahl's equation. 



<p(x)' 
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After the substitution y(x) = fi(r)(x) — 1) and the homogenizing substitution 
c) = ^gp, the Riccal 
to the canonical forms 



y(x) = 9 g ^ , the Riccati equation and it homogeneous counterpart are brought 



y'(x) = -y 2 (x) + y(x) ( - M ) + -y^ & g"(x) - z (x)g'(x) -u{x)g{x) = 0, 

(19) 

where 

Note that when q = 0, equation (fT9)l becomes essentially of first order g"(x) — 
g'(x)z(x) = 0, and <?'(x) = e z ^ x \ with Z(x) = J z{x)dx, recovering Segerdahl's 
result, see Example [21 

Note that having an explicit general solution 77(0;) to the Riccati equation of 
the system (|18l) leads to an explicit general solution of the system, obtained by 



M(x) = iexp ^ 



■q{t)dt — fj,x ] , 



where L is an arbitrary constant. Thus, the solution of the first-passage problem 
will be available analytically (up to quadratures), whenever the Riccati solution 
is. In other words, we have reduced the problem of solving system (fT2|) to solve 
a Riccati equation. Additionally, it is worth to remark that, in a similar way, 
solving higher order dimensional systems © can be also reduced to solving 
matrix Riccati equations. 

Riccati equations are Lie systems generally related to a Vessiot-Guldberg 
Lie algebra isomorphic to st(2,R), see [5]. It follows that this Lie algebra is not 
solvable, and the method exposed in [J does not apply. Nevertheless, the theory 
of Lie systems provides the following result allowing us to integrate Riccati 
equations for certain specific cases of their x-dependent coefficients. This result 
motivates the use of the alternative approach to systems (|T2l carried out in this 
Section. Although a similar result could also be demonstrated for these systems, 
it much simpler to use an alternative approach and use the below result, which 
proof can be found in [TO] . 

Theorem 1. The necessary and sufficient condition for the existence of a trans- 
formation 

f, = G(x)r], G(x) > 0, 

relating the Riccati equation 

P- = b (x) + bt(x)ri + b 2 (x)ri 2 , b b 2 ± 0, (21) 
dx 

to an integrable one given by 
dfj 



D(x)(c + Cl fj + c 2 f) , c c 2 =^0, (22) 
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where cq,ci,c 2 are real numbers and D(x) is a non-vanishing function, are 

where k = sg(D) = sgipa/co). The transformation is then uniquely defined by 

b 2 (x)c 
V bo(x)c 2 

As Riccati equations of the form (|22|) can be transformed easily into an 
autonomous equation by means of a x-reparametrization and these equations 
are always integrable by one quadrature, it follows that all the Riccati equations 
satisfying the conditions ([23)) can be solved. The family of Riccati equations 
satisfying the previous theorem slightly generalises some previous results on the 
topic. Indeed, the similarity with the result given by Allen and Stein in [2] is 
what motivated calling this family the generalised Allen-Stein family. 

Our main aim now is to determine which Segerdahl's equations can be in- 
tegrated by means of the above proposition in order to provide new models 
admitting exact solutions. In order to do so, note that Segerdahl's equations can 
be cast into the form (|2"Tj) . with 

&o(a0 = — h(x)=(n+^), h(x) = -n. (24) 

Substituting the above functions in the integrability condition (|23p , we get that 
Riccati equation (|18p is integrable if there exist constants Cq , C\ , and c 2 such that 
the drift ip(x) satisfies the equation 

(f'/2 + (X + q) + fMf = KCiy/ — fl\coC2ip. (25) 

For example, in the particular case C\ = 0, the above integrability condition 
reads 

ip' + 2fnp + 2 (A + q) = 0, (26) 
whose general solution, ipj((x), is 

ipK ( x ) = ±±l(Ke-^-l), (27) 
H 

with K a nonzero real constant. An explicit solution for the classical ruin problem 
with this drift follows. 



5 The Ruin probability 

We consider now the ruin problem with downwards exponential jumps and drift 
(|27[) . In the absence of jumps, we would have a dynamical system an attractive 

• , ln(K) 

point x = 2^T- 
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The most interesting exit problem is finding 

&(x) = P x [T L <n,T L <e q ], (28) 

where e q is an independent exponential horizon, i.e. P[e q > x] = e~ qx , and 
where I < L < x$. Note that in this case the jumps and drift go in opposite 
directions and both boundaries may be attained; therefore, this requires solving 
the system (|12[) with the boundary conditions M(l) = 0,!Z'(L) = 1. 

We turn now to the easiest case xo < I < with the drift negative, in which 
case the problem ([2"5]) (with M{1) = 0,!f r (L) = 1) has no solution, since because 
of the negative drift, with 'no vehicle going up', exit above is impossible. 

The other two sided problem 

&(x) = P x [T l <T L ,T l <e q \, (29) 

with boundary conditions M(0) = S'(O) = 1, coincides now with a one sided 
exit problem (again, since hitting the upper boundary is impossible). Sine the 
upper boundary is irrelevant, we will assume thus w.l.o.g. that L = oo and / = 0, 
or, equivalently, that K < 1 (note that <pk(0) = (K — 1) , lim-r^oo ipk{x) = 

M 

Under the above assumptions, the corresponding Segerdahl's equation be- 
comes integrable (by quadratures) by applying Theorem [T] with cq = 1, C2 = — 1, 
and ci = 0. More specifically, the change of variables given by 



V=J ^rfj, (30) 



transforms the Riccati equation in (|18[) into 

■(i-/r). (3D 



dx y <pk(x) 

In order to transform the Riccati equation (|31[) into an autonomous one, we 
define the reparametrisation 



-A/i _, , 1 / A , ( 1 - VT^K VI - e-^K + 1 
ax = 1,1 7~t^ x = ^ x ( x ) = ~\ r "T 



<Pk{x) 2\jq + \ - K+ 1 1 - VI - e- 2 ^Ji" 

In consequence, x(0) = and limx-j-oo x(x) = oo. Now, the solution for equation 
(|3Tjt reads 

e 2s -J^ 

with K\ a constant. Hence, in view of the change of variables (IBT)]) . we obtain that 
the general solution of the Riccati equation in (|18p under our present assumptions 
is 
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For each value of the function r)(x), we obtain, in view of the equation dM/dx 
(firj(x) — fi)M, that the function M{x) reads 

M{x) = Lexp (fx J n(x)dx - fix) = ^e'" x (K ie ~ x( - X ^ + e x{x) ). 

On the other hand, \P{x) = rj(x)M(x). Therefore, we obtain 



Summing up, the solutions for the system (|18[) in the particular case admitting 
the drift (J57J) are 



2 y q + A 
M{x) = -ye^^A-je^W + e x{x) ). 
Assume now that <f(0) = M(0) = 1. Therefore, 



(l + ifi(-ff)) y ? + A ^ 
where 



In order to analyse the behaviour of the above solutions at infinity, it is necessary 
to obtain an approximate expression of x[x) at x » 0. As y/1 + y ~ 1 + y/2 
when y ~ 0, then 



X W - W 77 — TT7 lo g 



4(g + A) b \ |1- Vl-e- 2 ^^! 



4(g + A) b Ve-^l^iy V A + 9 



Using the above result, it follows that 
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Consequently, it turns out that lim^^oo ty(x) = lim^^oo M(x) = 0. 

In order to illustrate the behaviour of the ruin probabilities ()32j) . we provide 
the following figure showing a particular case of our solutions and its correspond- 
ing drift lpk(x). 




4 6 8 10 



Fig. 1. Ruin probabilities and drift for /i = 1.5, A = q = 1/2, and K — 0.75. 
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